A Genomic Analysis of Esophageal Squamous Cell Carcinoma in Eastern Africa.

BACKGROUND: Esophageal squamous cell carcinoma (ESCC) comprises 90% of all esophageal cancer cases globally and is the most common histology in low-resource settings. Eastern Africa has a disproportionately high incidence of ESCC. METHODS: We describe the genomic profiles of 61 ESCC cases from Tanzania and compare them to profiles from an existing cohort of ESCC cases from Malawi. We also provide a comparison to ESCC tumors in The Cancer Genome Atlas (TCGA). RESULTS: We observed substantial transcriptional overlap with other squamous histologies via comparison with TCGA PanCan dataset. DNA analysis revealed known mutational patterns, both genome-wide as well as in genes known to be commonly mutated in ESCC. TP53 mutations were the most common somatic mutation in tumors from both Tanzania and Malawi but were detected at lower frequencies than previously reported in ESCC cases from other settings. In a combined analysis, two unique transcriptional clusters were identified: a proliferative/epithelial cluster and an invasive/migrative/mesenchymal cluster. Mutational signature analysis of the Tanzanian cohort revealed common signatures associated with aging and cytidine deaminase activity (APOBEC) and an absence of signature 29, which was previously reported in the Malawi cohort. CONCLUSIONS: This study defines the molecular characteristics of ESCC in Tanzania, and enriches the Eastern African dataset, with findings of overall similarities but also some heterogeneity across two unique sites. IMPACT: Despite a high burden of ESCC in Eastern Africa, investigations into the genomics in this region are nascent. This represents the largest comprehensive genomic analysis ESCC from sub-Saharan Africa to date.

Inferred Immune-Cell Activity Is an Independent Predictor of HER2-Negative Breast Cancer Prognosis and Response to Paclitaxel-Based Therapy in the GeparSepto Trial.

PURPOSE: Tumor microenvironment (TME) immune markers have been correlated with both response to neoadjuvant therapy and prognosis in patients with breast cancer. Here, immune-cell activity of breast cancer tumors was inferred by expression-based analysis to determine if it is prognostic and/or predictive of response to neoadjuvant paclitaxel-based therapy in the GeparSepto (G7) trial (NCT01583426). EXPERIMENTAL DESIGN: Pre-study biopsies from 279 patients with HER2-negative breast cancer in the G7 trial underwent RNA-seq-based profiling of 104 immune-cell-specific genes to assess inferred Immune Cell Activity (iICA) of 23 immune-cell types. Hierarchical clustering was used to classify tumors as iICA "hot," "warm," or "cold" by comparison of iICA in the G7 cohort relative to that of 1,467 samples from a tumor database established by Nantomics LLC. Correlations between iICA cluster, pathology-assessed tumor-infiltrating lymphocytes (TIL), and hormone receptor (HR) status for pathologic complete response (pCR), disease-free survival (DFS), and overall survival (OS) were determined. RESULTS: iICA cluster correlated with TIL levels. The highest pCR rates were observed in hot cluster tumors, and those with relatively higher TILs. Greater inferred activity of several T-cell types was significantly associated with pCR and survival. DFS and OS were prolonged in patients with hot or warm cluster tumors, the latter particularly for HR negative tumors, even if TILs were relatively low. CONCLUSIONS: Overall, TIL level better predicted pCR, but iICA cluster better predicted survival. Differences in associations between TILs, cluster, pCR, and survival were observed for HR-positive tumors versus HR-negative tumors, suggesting expanded study of the implication of these findings is warranted.

Case report: Complete pathologic response with first-line immunotherapy combination in a young adult with massive liver dissemination of mismatch repair-deficient metastatic colorectal cancer: Immunological and molecular profiling.

The current level of evidence for immunotherapy in previously untreated microsatellite unstable metastatic colorectal cancer is based on recent pieces of evidence of few studies that demonstrated durable response and clinical benefit, in terms of objective response rate, disease control rate, and progression-free survival in this subgroup of patients. On the basis of combinatorial immunotherapy with nivolumab plus ipilimumab, we report the exceptional case of a complete pathological response in a 21-year-old woman presenting a clinically aggressive stage IV colorectal cancer with massive nodal and liver involvement. Extensive molecular analyses based on whole genome next-generation DNA sequencing, RNA sequencing, fluorescent multiplex immunohistochemistry, and flow cytometry provided a detailed description of tumoral and immunological characteristics of this noteworthy clinical case.

A single-cell transcriptional gradient in human cutaneous memory T cells restricts Th17/Tc17 identity.

The homeostatic mechanisms that fail to restrain chronic tissue inflammation in diseases, such as psoriasis vulgaris, remain incompletely understood. We profiled transcriptomes and epitopes of single psoriatic and normal skin-resident T cells, revealing a gradated transcriptional program of coordinately regulated inflammation-suppressive genes. This program, which is sharply suppressed in lesional skin, strikingly restricts Th17/Tc17 cytokine and other inflammatory mediators on the single-cell level. CRISPR-based deactivation of two core components of this inflammation-suppressive program, ZFP36L2 and ZFP36, replicates the interleukin-17A (IL-17A), granulocyte macrophage-colony-stimulating factor (GM-CSF), and interferon gamma (IFNγ) elevation in psoriatic memory T cells deficient in these transcripts, functionally validating their influence. Combinatoric expression analysis indicates the suppression of specific inflammatory mediators by individual program members. Finally, we find that therapeutic IL-23 blockade reduces Th17/Tc17 cell frequency in lesional skin but fails to normalize this inflammatory-suppressive program, suggesting how treated lesions may be primed for recurrence after withdrawal of treatment.

Safety, Feasibility, and Merits of Longitudinal Molecular Testing of Multiple Metastatic Sites to Inform mTNBC Patient Treatment in the Intensive Trial of Omics in Cancer.

PURPOSE: Patients with metastatic triple-negative breast cancer (mTNBC) have poor outcomes. The Intensive Trial of Omics in Cancer (ITOMIC) sought to determine the feasibility and potential efficacy of informing treatment decisions through multiple biopsies of mTNBC deposits longitudinally over time, accompanied by analysis using a distributed network of experts. METHODS: Thirty-one subjects were enrolled and 432 postenrollment biopsies performed (clinical and study-directed) of which 332 were study-directed. Molecular profiling included whole-genome sequencing or whole-exome sequencing, cancer-associated gene panel sequencing, RNA-sequencing, and immunohistochemistry. To afford time for analysis, subjects were initially treated with cisplatin (19 subjects), or another treatment they had not received previously. The results were discussed at a multi-institutional ITOMIC Tumor Board, and a report transmitted to the subject's oncologist who arrived at the final treatment decision in conjunction with the subject. Assistance was provided to access treatments that were predicted to be effective. RESULTS: Multiple biopsies in single settings and over time were safe, and comprehensive analysis was feasible. Two subjects were found to have lung cancer, one had carcinoma of unknown primary site, tumor samples from three subjects were estrogen receptor-positive and from two others, human epidermal growth factor receptor 2-positive. Two subjects withdrew. Thirty-four of 112 recommended treatments were accessed using approved drugs, clinical trials, and single-patient investigational new drugs. After excluding the three subjects with nonbreast cancers and the two subjects who withdrew, 22 of 26 subjects (84.6%) received at least one ITOMIC Tumor Board-recommended treatment. CONCLUSION: Further exploration of this approach in patients with mTNBC is merited.

Identification and validation of expressed HLA-binding breast cancer neoepitopes for potential use in individualized cancer therapy.

BACKGROUND: Therapeutic regimens designed to augment the immunological response of a patient with breast cancer (BC) to tumor tissue are critically informed by tumor mutational burden and the antigenicity of expressed neoepitopes. Herein we describe a neoepitope and cognate neoepitope-reactive T-cell identification and validation program that supports the development of next-generation immunotherapies. METHODS: Using GPS Cancer, NantOmics research, and The Cancer Genome Atlas databases, we developed a novel bioinformatic-based approach which assesses mutational load, neoepitope expression, human leukocyte antigen (HLA)-binding prediction, and in vitro confirmation of T-cell recognition to preferentially identify targetable neoepitopes. This program was validated by application to a BC cell line and confirmed using tumor biopsies from two patients with BC enrolled in the Tumor-Infiltrating Lymphocytes and Genomics (TILGen) study. RESULTS: The antigenicity and HLA-A2 restriction of the BC cell line predicted neoepitopes were determined by reactivity of T cells from HLA-A2-expressing healthy donors. For the TILGen subjects, tumor-infiltrating lymphocytes (TILs) recognized the predicted neoepitopes both as peptides and on retroviral expression in HLA-matched Epstein-Barr virus-lymphoblastoid cell line and BC cell line MCF-7 cells; PCR clonotyping revealed the presence of T cells in the periphery with T-cell receptors for the predicted neoepitopes. These high-avidity immune responses were polyclonal, mutation-specific and restricted to either HLA class I or II. Interestingly, we observed the persistence and expansion of polyclonal T-cell responses following neoadjuvant chemotherapy. CONCLUSIONS: We demonstrate our neoepitope prediction program allows for the successful identification of neoepitopes targeted by TILs in patients with BC, providing a means to identify tumor-specific immunogenic targets for individualized treatment, including vaccines or adoptively transferred cellular therapies.

Single-Cell Profiling Reveals Divergent, Globally Patterned Immune Responses in Murine Skin Inflammation.

Inflammatory response heterogeneity has impeded high-resolution dissection of diverse immune cell populations during activation. We characterize mouse cutaneous immune cells by single-cell RNA sequencing, after inducing inflammation using imiquimod and oxazolone dermatitis models. We identify 13 CD45+ subpopulations, which broadly represent most functionally characterized immune cell types. Oxazolone pervasively upregulates Jak2/Stat3 expression across T cells and antigen-presenting cells (APCs). Oxazolone also induces Il4/Il13 expression in newly infiltrating basophils, and Il4ra and Ccl24, most prominently in APCs. In contrast, imiquimod broadly upregulates Il17/Il22 and Ccl4/Ccl5. A comparative analysis of single-cell inflammatory transcriptional responses reveals that APC response to oxazolone is tightly restricted by cell identity, whereas imiquimod enforces shared programs on multiple APC populations in parallel. These global molecular patterns not only contrast immune responses on a systems level but also suggest that the mechanisms of new sources of inflammation can eventually be deduced by comparison to known signatures.

Large scale, robust, and accurate whole transcriptome profiling from clinical formalin-fixed paraffin-embedded samples.

Transcriptome profiling can provide information of great value in clinical decision-making, yet RNA from readily available formalin-fixed paraffin-embedded (FFPE) tissue is often too degraded for quality sequencing. To assess the clinical utility of FFPE-derived RNA, we performed ribo-deplete RNA extractions on > 3200 FFPE slide samples; 25 of these had direct FFPE vs. fresh frozen (FF) replicates, 57 were sequenced in 2 different labs, 87 underwent multiple library analyses, and 16 had direct microdissected vs. macrodissected replicates. Poly-A versus ribo-depletion RNA extraction methods were compared using transcriptomes of TCGA cohort and 3116 FFPE samples. Compared to FF, FFPE transcripts coding for nuclear/cytoplasmic proteins involved in DNA packaging, replication, and protein synthesis were detected at lower rates and zinc finger family transcripts were of poorer quality. The greatest difference in extraction methods was in histone transcripts which typically lack poly-A tails. Encouragingly, the overall sequencing success rate was 81%. Exome coverage was highly concordant in direct FFPE and FF replicates, with 98% agreement in coding exon coverage and a median correlation of whole transcriptome profiles of 0.95. We provide strong rationale for clinical use of FFPE-derived RNA based on the robustness, reproducibility, and consistency of whole transcriptome profiling.

A risk-associated Active transcriptome phenotype expressed by histologically normal human breast tissue and linked to a pro-tumorigenic adipocyte population.

BACKGROUND: Previous studies have identified and validated a risk-associated Active transcriptome phenotype commonly expressed in the cancer-adjacent and histologically normal epithelium, stroma, and adipose containing peritumor microenvironment of clinically established invasive breast cancers, conferring a 2.5- to 3-fold later risk of dying from recurrent breast cancer. Expression of this Active transcriptome phenotype has not yet been evaluated in normal breast tissue samples unassociated with any benign or malignant lesions; however, it has been associated with increased peritumor adipocyte composition. METHODS: Detailed histologic and transcriptomic (RNAseq) analyses were performed on normal breast biopsy samples from 151 healthy, parous, non-obese (mean BMI = 29.60 ± 7.92) women, ages 27-66 who donated core breast biopsy samples to the Komen Tissue Bank, and whose average breast cancer risk estimate (Gail score) at the time of biopsy (1.27 ± 1.34) would not qualify them for endocrine prevention therapy. RESULTS: Full genome RNA sequencing (RNAseq) identified 52% (78/151) of these normal breast samples as expressing the Active breast phenotype. While Active signature genes were found to be most variably expressed in mammary adipocytes, donors with the Active phenotype had no difference in BMI but significantly higher Gail scores (1.46 vs. 1.18; p = 0.007). Active breast samples possessed 1.6-fold more (~ 80%) adipocyte nuclei, larger cross-sectional adipocyte areas (p < 0.01), and 0.5-fold fewer stromal and epithelial cell nuclei (p < 1e-6). Infrequent low-level expression of cancer gene hotspot mutations was detected but not enriched in the Active breast samples. Active samples were enriched in gene sets associated with adipogenesis and fat metabolism (FDR q ≤ 10%), higher signature scores for cAMP-dependent lipolysis known to drive breast cancer progression, white adipose tissue browning (Wilcoxon p < 0.01), and genes associated with adipocyte activation (leptin, adiponectin) and remodeling (CAV1, BNIP3), adipokine growth factors (IGF-1, FGF2), and pro-inflammatory fat signaling (IKBKG, CCL13). CONCLUSIONS: The risk-associated Active transcriptome phenotype first identified in cancer-adjacent breast tissues also occurs commonly in healthy women without breast disease who do not qualify for breast cancer chemoprevention, and independently of breast expressed cancer-associated mutations. The risk-associated Active phenotype appears driven by a pro-tumorigenic adipocyte microenvironment that can predate breast cancer development.

Transcriptomic silencing as a potential mechanism of treatment resistance.

Next-generation sequencing (NGS) has not revealed all the mechanisms underlying resistance to genomically matched drugs. Here, we performed in 1417 tumors whole-exome tumor (somatic)/normal (germline) NGS and whole-transcriptome sequencing, the latter focusing on a clinically oriented 50-gene panel in order to examine transcriptomic silencing of putative driver alterations. In this large-scale study, approximately 13% of the somatic single nucleotide variants (SNVs) were unexpectedly not expressed as RNA; 23% of patients had ≥1 nonexpressed SNV. SNV-bearing genes consistently transcribed were TP53, PIK3CA, and KRAS; those with lower transcription rates were ALK, CSF1R, ERBB4, FLT3, GNAS, HNF1A, KDR, PDGFRA, RET, and SMO. We also determined the frequency of tumor mutations being germline, rather than somatic, in these and an additional 462 tumors with tumor/normal exomes; 33.8% of germline SNVs within the gene panel were rare (not found after filtering through variant information domains) and at risk of being falsely reported as somatic. Both the frequency of silenced variant transcription and the risk of falsely identifying germline mutations as somatic/tumor related are important phenomena. Therefore, transcriptomics is a critical adjunct to genomics when interrogating patient tumors for actionable alterations, because, without expression of the target aberrations, there will likely be therapeutic resistance.

Overcoming hypoxia-induced functional suppression of NK cells.

BACKGROUND: Natural killer (NK) cells are immune cells capable of killing virally infected cells and tumor cells without the need for antigen stimulation. Tumors, however, can create a suppressive microenvironment that decreases NK function. A feature of many tumors is hypoxia (low oxygen perfusion), which has been previously shown to decrease NK function. A high affinity NK (haNK) cell has been engineered to express a high affinity CD16 receptor as well as internal interleukin (IL)-2 for increased antibody-dependent cellular cytotoxicity (ADCC) and activation, respectively. We sought to investigate the tolerance of NK cells versus haNK cells to hypoxia. METHODS: We exposed healthy donor (HD) NK and X-irradiated haNK cells to normoxia (20% oxygen) as well as hypoxia (0% oxygen) and investigated their ability to kill prostate, breast and lung tumor cell lines after 5 hours. We also used monoclonal antibodies cetuximab (anti-EGFR) or avelumab (antiprogrammed death-ligand 1) to investigate the effects of hypoxia on NK ADCC. Genomic and proteomic analyzes were done to determine the effect of hypoxia on the expression of factors important to NK cell function. RESULTS: While HD NK cell cytolytic abilities were markedly and significantly impaired under hypoxic conditions, haNK cells maintained killing capacity under hypoxic conditions. NK killing, serial killing and ADCC were maintained under hypoxia in haNK cells. IL-2 has been previously implicated in serial killing and perforin regeneration and thus the endogenous IL-2 produced by haNK cells is likely a driver of the maintained killing capacity of haNK cells under hypoxic conditions. Activation of signal transducer and activator of transcription 3 (STAT3) is not seen in haNKs under hypoxia but is significant in HD NK cells. Pharmaceutical activation of STAT3 in haNKs led to reduced killing, implicating active STAT3 in reduced NK cell function. CONCLUSIONS: In contrast to HD NK cells, haNK cells are resistant to acute hypoxia. The potent cytolytic function of haNK cells was maintained in an environment comparable to what would be encountered in a tumor. The data presented here provide an additional mechanism of action for haNK cells that are currently being evaluated in clinical trials for several tumor types.

A deep learning image-based intrinsic molecular subtype classifier of breast tumors reveals tumor heterogeneity that may affect survival.

BACKGROUND: Breast cancer intrinsic molecular subtype (IMS) as classified by the expression-based PAM50 assay is considered a strong prognostic feature, even when controlled for by standard clinicopathological features such as age, grade, and nodal status, yet the molecular testing required to elucidate these subtypes is not routinely performed. Furthermore, when such bulk assays as RNA sequencing are performed, intratumoral heterogeneity that may affect prognosis and therapeutic decision-making can be missed. METHODS: As a more facile and readily available method for determining IMS in breast cancer, we developed a deep learning approach for approximating PAM50 intrinsic subtyping using only whole-slide images of H&E-stained breast biopsy tissue sections. This algorithm was trained on images from 443 tumors that had previously undergone PAM50 subtyping to classify small patches of the images into four major molecular subtypes-Basal-like, HER2-enriched, Luminal A, and Luminal B-as well as Basal vs. non-Basal. The algorithm was subsequently used for subtype classification of a held-out set of 222 tumors. RESULTS: This deep learning image-based classifier correctly subtyped the majority of samples in the held-out set of tumors. However, in many cases, significant heterogeneity was observed in assigned subtypes across patches from within a single whole-slide image. We performed further analysis of heterogeneity, focusing on contrasting Luminal A and Basal-like subtypes because classifications from our deep learning algorithm-similar to PAM50-are associated with significant differences in survival between these two subtypes. Patients with tumors classified as heterogeneous were found to have survival intermediate between Luminal A and Basal patients, as well as more varied levels of hormone receptor expression patterns. CONCLUSIONS: Here, we present a method for minimizing manual work required to identify cancer-rich patches among all multiscale patches in H&E-stained WSIs that can be generalized to any indication. These results suggest that advanced deep machine learning methods that use only routinely collected whole-slide images can approximate RNA-seq-based molecular tests such as PAM50 and, importantly, may increase detection of heterogeneous tumors that may require more detailed subtype analysis.

Multiregion exome sequencing of ovarian immature teratomas reveals 2N near-diploid genomes, paucity of somatic mutations, and extensive allelic imbalances shared across mature, immature, and disseminated components.

Immature teratoma is a subtype of malignant germ cell tumor of the ovary that occurs most commonly in the first three decades of life, frequently with bilateral ovarian disease. Despite being the second most common malignant germ cell tumor of the ovary, little is known about its genetic underpinnings. Here we performed multiregion whole-exome sequencing to interrogate the genetic zygosity, clonal relationship, DNA copy number, and mutational status of 52 pathologically distinct tumor components from ten females with ovarian immature teratomas, with bilateral tumors present in five cases and peritoneal dissemination in seven cases. We found that ovarian immature teratomas are genetically characterized by 2N near-diploid genomes with extensive loss of heterozygosity and an absence of genes harboring recurrent somatic mutations or known oncogenic variants. All components within a single ovarian tumor (immature teratoma, mature teratoma with different histologic patterns of differentiation, and yolk sac tumor) were found to harbor an identical pattern of loss of heterozygosity across the genome, indicating a shared clonal origin. In contrast, the four analyzed bilateral teratomas showed distinct patterns of zygosity changes in the right versus left sided tumors, indicating independent clonal origins. All disseminated teratoma components within the peritoneum (including gliomatosis peritonei) shared a clonal pattern of loss of heterozygosity with either the right or left primary ovarian tumor. The observed genomic loss of heterozygosity patterns indicate that diverse meiotic errors contribute to the formation of ovarian immature teratomas, with 11 out of the 15 genetically distinct clones determined to result from nondisjunction errors during meiosis I or II. Overall, these findings suggest that copy-neutral loss of heterozygosity resulting from meiotic abnormalities may be sufficient to generate ovarian immature teratomas from germ cells.

Efficient Tumor Clearance and Diversified Immunity through Neoepitope Vaccines and Combinatorial Immunotherapy.

Progressive tumor growth is associated with deficits in the immunity generated against tumor antigens. Vaccines targeting tumor neoepitopes have the potential to address qualitative defects; however, additional mechanisms of immune failure may underlie tumor progression. In such cases, patients would benefit from additional immune-oncology agents targeting potential mechanisms of immune failure. This study explores the identification of neoepitopes in the MC38 colon carcinoma model by comparison of tumor to normal DNA and tumor RNA sequencing technology, as well as neoepitope delivery by both peptide- and adenovirus-based vaccination strategies. To improve antitumor efficacies, we combined the vaccine with a group of rationally selected immune-oncology agents. We utilized an IL15 superagonist to enhance the development of antigen-specific immunity initiated by the neoepitope vaccine, PD-L1 blockade to reduce tumor immunosuppression, and a tumor-targeted IL12 molecule to facilitate T-cell function within the tumor microenvironment. Analysis of tumor-infiltrating leukocytes demonstrated this multifaceted treatment regimen was required to promote the influx of CD8+ T cells and enhance the expression of transcripts relating to T-cell activation/effector function. Tumor-targeted IL12 resulted in a marked increase in clonality of T-cell repertoire infiltrating the tumor, which when sculpted with the addition of either a peptide or adenoviral neoepitope vaccine promoted efficient tumor clearance. In addition, the neoepitope vaccine induced the spread of immunity to neoepitopes expressed by the tumor but not contained within the vaccine. These results demonstrate the importance of combining neoepitope-targeting vaccines with a multifaceted treatment regimen to generate effective antitumor immunity.

Loss of ZNF750 in ocular and cutaneous sebaceous carcinoma.

BACKGROUND: Sebaceous carcinoma (SeC) is an uncommon malignancy arising from sebaceous glands of the conjunctiva and skin. Recurrent mutations in the ZNF750 were recently identified in ocular SeC. We assessed whether ZNF750 loss is a specific feature of ocular SeC or a general feature of sebaceous tumors. METHODS: Immunostaining for ZNF750 expression was performed in 54 benign and malignant sebocytic proliferations. Staining for ZNF750 was scored on a three-tier scale: positive (>75%), partially positive (5%-74%), and negative (<5%). RESULTS: ZNF750 expression was negative in 4/11 ocular SeC, and partially positive in 4/11 ocular SeC and 6/13 cutaneous SeC. No extraocular tumors were negative. No loss was found in sebaceous adenoma or sebaceous hyperplasia. In nine previously sequenced ocular SeCs, two lacked detectable somatic mutations in ZNF750, but showed complete loss of staining, indicating non-mutational inactivation of ZNF750. CONCLUSION: We show complete loss of the ZNF750 epidermal differentiation regulator in about half of ocular SeC, highlighting the most common genetic defect in this cancer type. Loss of ZNF750 expression is seen even in tumors without truncating mutations and reduced in many of the remaining ocular and cutaneous SeC. In contrast, no ZNF750 loss was detected in benign sebaceous proliferations.

A20 and ABIN1 Suppression of a Keratinocyte Inflammatory Program with a Shared Single-Cell Expression Signature in Diverse Human Rashes.

Genetic variation in the NF-κB inhibitors, ABIN1 and A20, increase risk for psoriasis. While critical for hematopoietic immune cell function, these genes are believed to additionally inhibit psoriasis by dampening inflammatory signaling in keratinocytes. We dissected ABIN1 and A20's regulatory role in human keratinocyte inflammation using an RNA sequencing-based comparative genomic approach. Here we show subsets of the IL-17 and tumor necrosis factor-α signaling pathways are robustly restricted by A20 overexpression. In contrast, ABIN1 overexpression inhibits these genes more modestly for IL-17, and weakly for tumor necrosis factor-α. Our genome-scale analysis also indicates that inflammatory program suppression appears to be the major transcriptional influence of A20/ABIN1 overexpression, without obvious influence on keratinocyte viability genes. Our findings thus enable dissection of the differing anti-inflammatory mechanisms of two distinct psoriasis modifiers, which may be directly exploited for therapeutic purposes. Importantly, we report that IL-17-induced targets of A20 show similar aberrant epidermal layer-specific transcriptional upregulation in keratinocytes from diseases as diverse as psoriasis, atopic dermatitis, and erythrokeratodermia variabilis, suggesting a contributory role for epidermal inflammation in a broad spectrum of rashes.

Transcriptional Programming of Normal and Inflamed Human Epidermis at Single-Cell Resolution.

Perturbations in the transcriptional programs specifying epidermal differentiation cause diverse skin pathologies ranging from impaired barrier function to inflammatory skin disease. However, the global scope and organization of this complex cellular program remain undefined. Here we report single-cell RNA sequencing profiles of 92,889 human epidermal cells from 9 normal and 3 inflamed skin samples. Transcriptomics-derived keratinocyte subpopulations reflect classic epidermal strata but also sharply compartmentalize epithelial functions such as cell-cell communication, inflammation, and WNT pathway modulation. In keratinocytes, ∼12% of assessed transcript expression varies in coordinate patterns, revealing undescribed gene expression programs governing epidermal homeostasis. We also identify molecular fingerprints of inflammatory skin states, including S100 activation in the interfollicular epidermis of normal scalp, enrichment of a CD1C+CD301A+ myeloid dendritic cell population in psoriatic epidermis, and IL1ßhiCCL3hiCD14+ monocyte-derived macrophages enriched in foreskin. This compendium of RNA profiles provides a critical step toward elucidating epidermal diseases of development, differentiation, and inflammation.

APOBEC mutation drives early-onset squamous cell carcinomas in recessive dystrophic epidermolysis bullosa.

Recessive dystrophic epidermolysis bullosa (RDEB) is a rare inherited skin and mucous membrane fragility disorder complicated by early-onset, highly malignant cutaneous squamous cell carcinomas (SCCs). The molecular etiology of RDEB SCC, which arises at sites of sustained tissue damage, is unknown. We performed detailed molecular analysis using whole-exome, whole-genome, and RNA sequencing of 27 RDEB SCC tumors, including multiple tumors from the same patient and multiple regions from five individual tumors. We report that driver mutations were shared with spontaneous, ultraviolet (UV) light-induced cutaneous SCC (UV SCC) and head and neck SCC (HNSCC) and did not explain the early presentation or aggressive nature of RDEB SCC. Instead, endogenous mutation processes associated with apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like (APOBEC) deaminases dominated RDEB SCC. APOBEC mutation signatures were enhanced throughout RDEB SCC tumor evolution, relative to spontaneous UV SCC and HNSCC mutation profiles. Sixty-seven percent of RDEB SCC driver mutations was found to emerge as a result of APOBEC and other endogenous mutational processes previously associated with age, potentially explaining a >1000-fold increased incidence and the early onset of these SCCs. Human papillomavirus-negative basal and mesenchymal subtypes of HNSCC harbored enhanced APOBEC mutational signatures and transcriptomes similar to those of RDEB SCC, suggesting that APOBEC deaminases drive other subtypes of SCC. Collectively, these data establish specific mutagenic mechanisms associated with chronic tissue damage. Our findings reveal a cause for cancers arising at sites of persistent inflammation and identify potential therapeutic avenues to treat RDEB SCC.

Identification of an immune gene expression signature associated with favorable clinical features in Treg-enriched patient tumor samples.

Immune heterogeneity within the tumor microenvironment undoubtedly adds several layers of complexity to our understanding of drug sensitivity and patient prognosis across various cancer types. Within the tumor microenvironment, immunogenicity is a favorable clinical feature in part driven by the antitumor activity of CD8+ T cells. However, tumors often inhibit this antitumor activity by exploiting the suppressive function of regulatory T cells (Tregs), thus suppressing the adaptive immune response. Despite the seemingly intuitive immunosuppressive biology of Tregs, prognostic studies have produced contradictory results regarding the relationship between Treg enrichment and survival. We therefore analyzed RNA-seq data of Treg-enriched tumor samples to derive a pan-cancer gene signature able to help reconcile the inconsistent results of Treg studies, by better understanding the variable clinical association of Tregs across alternative tumor contexts. We show that increased expression of a 32-gene signature in Treg-enriched tumor samples (n = 135) is able to distinguish a cohort of patients associated with chemosensitivity and overall survival. This cohort is also enriched for CD8+ T cell abundance, as well as the antitumor M1 macrophage subtype. With a subsequent validation in a larger TCGA pool of Treg-enriched patients (n = 626), our results reveal a gene signature able to produce unsupervised clusters of Treg-enriched patients, with one cluster of patients uniquely representative of an immunogenic tumor microenvironment. Ultimately, these results support the proposed gene signature as a putative biomarker to identify certain Treg-enriched patients with immunogenic tumors that are more likely to be associated with features of favorable clinical outcome.

Comprehensive genomic transcriptomic tumor-normal gene panel analysis for enhanced precision in patients with lung cancer.

A CMS approved test for lung cancer is based on tumor-only analysis of a targeted 35 gene panel, specifically excluding the use of the patient's normal germline tissue. However, this tumor-only approach increases the risk of mistakenly identifying germline single nucleotide polymorphisms (SNPs) as somatically-derived cancer driver mutations (false positives). 621 patients with 30 different cancer types, including lung cancer, were studied to compare the precision of tumor somatic variant calling in 35 genes using tumor-only DNA sequencing versus tumor-normal DNA plus RNA sequencing. When sequencing of lung cancer was performed using tumor genomes alone without normal germline controls, 94% of variants identified were SNPs and thus false positives. Filtering for common SNPs still resulted in as high as 48% false positive variant calling. With tumor-only sequencing, 29% of lung cancer patients had a false positive variant call in at least one of twelve genes with directly targetable drugs. RNA analysis showed 18% of true somatic variants were not expressed. Thus, sequencing and analysis of both normal germline and tumor genomes is necessary for accurate identification of molecular targets. Treatment decisions based on tumor-only analysis may result in the administration of ineffective therapies while also increasing the risk of negative drug-related side effects.